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Abstract. We discuss the general and approximate angular diameter distance in the Friedman-Robertson- Walker 
cosmological models with nonzero cosmological constant. We modify the equation for the angular diameter distance 
by taking into account the fact that locally the distribution of matter is non homogeneous. We present exact 
solutions of this equation in a few special cases. We propose an approximate analytic solution of this equation 
which is simple enough and sufficiently accurate to be useful in practical applications. 
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1. Introduction 

Recent observations of the type la supernovae and CMB 
anisotropy strongly indicate that the total matter-energy 
density of the universe is now dominated by some kind of 
vacuum energy also called " dark energy" or the cosmolog- 
ical constant A l|Perlmutter 19971 iPerlmutteret al.19991 
IRiess et al. 19981 IKiess 2000[) . The origin and nature 
of this vacuum energy remains unknown. There are 
several review articles providing thorough discussion 
of the history, interpretations and problems connected 
with the vacuum energy and observational constrains 
HZerdovich 19671 |Weinberg 1989|IHarroIl et a!.1992l> . 
The type la supernovae have been already observed at 
redshifts z > 1. It is well known from galaxy surveys 
that galaxies and clusters of galaxies up to a scale of 
~ lGpc are distributed non homogeneously forming 
filaments, walls and underdense voids. This indicates that 
on similar scales also the dark matter is distributed non 
homogeneously. In this paper we analyze the influence of 
local non homogeneities on the angular diameter distance 
in a universe with non zero cosmological constant. 
The angular diameter distance in a locally 
non homogeneous universe was discussed by 
Zeldovich ([Dashveski & Zcldovich 1965) see 
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also dDashveski fc Slysh 1966D and d Weinberg 1 989 
|Kayser et al. 1997). Later Dyer& Roeder (1972) used the 
so called empty beam approximation to derive an equa- 
tion for the angular diameter distance; for more detailed 
references see also HSchneider, Ehlers fc Falco~9 92), 
l|Kantowski 1998jl . and IjTomita fc Futamase 19991) . We 
follow the Dyer & Roeder method to derive the equation 
for angular diameter distance in a locally non homo- 
geneous universe with a cosmological constant. In the 
general case the equation for the angular diameter dis- 
tance does not have analytical solutions, it can be solved 
only numerically ( |Kantowski, Kao fc Thomas 2000| ). We 
have found an analytic approximate solution of this 
equation, which is simple and accurate enough to be 
useful in practical applications. This allowed us to find an 
approximate dependence of the angular diameter distance 
on the basic cosmological parameters. 
The paper is organized as follows: we begin with the 
general form of Sachs equations describing light prop- 
agation in an arbitrary spacetime and using the empty 
beam approximation we derive the equation for angular 
diameter distance. Then, we discuss the angular diameter 
distance as observed from a location not at the origin 
(z=0). Finally, after discussing some properties of the 
analytical solutions of the equation for the angular 
diameter distance in a locally non homogeneous universe, 
we propose an analytic approximation solution of this 
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equation valid in a wide redshift interval (0, 10). We 
have been motivated by the great advances in observing 
further and further objects and the construction of 
very powerful telescopes that provide the possibility 
to observe gravitational lensing by clusters of galaxies, 
and supernovae at large distances. It is therefore neces- 
sary to develop a more accurate formalism to describe 
the distance redshift relation and gravitational lensing 
by high-redshift objects. To achieve this goal we use 
cosmological models with realistic values of the basic 
cosmological parameters such as the Hubble constant and 
the average matter-energy density. In concluding remarks 
we summarize our results and discuss some perspectives . 

2. Propagation of light and local non 
homogeneities 

2.1. General considerations 

Let us consider a beam of light emanating from a source S 
in an arbitrary spacetime described by the metric tensor 
g a p. The light rays propagate along a null surface E, which 
is determined by the eikonal equation 



ff a/3 E, Q E )/3 = 0. 



(1) 



A light ray is identified with a null geodesic on E with 
the tangent vector k a = — E, Q . The light rays in the beam 
can be described by x a = x a (v,y a ), where v is an affine 
parameter, and y a (a=l, 2, 3) are three parameters spec- 
ifying different rays. The vector field tangent to the light 
dx a 

ray congruence, k a — — — = — E, a , determines two op- 
era 

tical scalars, the expansion 6 and the shear a, which are 
defined by 



-k a ; a , a = k a .f3m a m 13 , 



(2) 



where m a = -^(£ a — ii] a ) is a complex vector span- 
ning the spacelike 2-space (the screen space) orthogonal 
to k a (k a fh a = 0). Since k a = — E, Q the vorticity con- 
nected with the light beam is zero in all our considera- 
tions; therefore in this case a and 8 fully characterize the 
congruence. These two optical scalars satisfy the Sachs 
( Sachs & Kristian 1966) propagation equations 



9 2 + \a\ 2 = -lR a0 k a k 



(3) 



a + 26a = -^C aM 8m a k l3 m<k s , (4) 

where dot denotes the derivative with respect to v, R a p 
is the Ricci tensor, and C a0J s is the Weyl tensor. Eqs. 
JSJ and Q follow from the Ricci identity. The optical 
scalars and a describe the relative rate of change of an 
infinitesimal area A of the cross section of the beam of 
light rays and its distortion. The expansion 9 is related to 
the relative change of an infinitesimal surface area A of 
the beam's cross section by 



Let us use these equations to study the propagation of 
light in the Friedman-Robertson- Walker (FRW) space- 
time. The FRW spacetime is conformally flat, so that 
C a /3 7 s = 0. From Eq. |T3J| it follows that, in the FRW 
spacetime if the shear of the null ray congruence is ini- 
tially zero, then it always vanishes. Therefore, assuming 
that the light beam emanating from the source S has van- 
ishing shear, we can disregard the shear parameter alto- 
gether (empty beam approximation) . Using Eq. J^J we can 
rewrite equation J2J in the form 

VA+^R afi k a k p VA = 0. (6) 

An observer moving with the 4- velocity vector u a will as- 
sociate with the light-ray a circular frequency lo — cu a k a . 
Different observers assign different frequencies to the same 
light ray. The shift of frequency, as measured by an ob- 
server comoving with the source and an arbitrary observer, 
is related to their relative velocity or the redshift z by 



1 + z = — = —u a k° 



(7) 



where u> is the frequency measured by a comoving observer 
and wo by an arbitrary observer. Differentiating this equa- 
tion with respect to the affine parameter v, we obtain 

dz 



dv 



-k a k$ 



U a -/3. 



(8) 



Since the angular diameter distance D is proportional to 
y/A we can rewrite Eq.® using D instead of y/~A and, 
at the same time, replace the affine parameter v by the 
redshift z; we obtain 



dz\ 2 d 2 D 
dv J dz 2 



d 2 z\ dD 
dv 2 I dz 



A-kG 



+ -^-T a0 k a k p D = 0, (9) 



where we have used the Einstein equations to replace the 
Ricci tensor by the energy-momentum tensor. A solution 
of Eq. © is related to the angular diameter distance if it 
satisfies the following initial conditions: 



2?(*)U=o = 0, 



dD{z) | c 
dz ~ Hq 



(10) 



2.2. The angular diameter distance between two 
objects at different redshifts 

In order to use solutions of the Eq. 10 to describe gravi- 
tational lenses we have to introduce the angular diameter 
distance between the source and the lens D(zi, z s ), where 
zi and z s denote correspondingly the redshifts of the lens 
and the source. Let D(z\, z%) denote the angular diameter 
distance between a fictitious observer at z\ and a source 
at z 2 ; of course, -D(0, z) — D(z). Suppose that we know 
the general solution of Eq. @ for D(z) 1 which satisfies the 
initial conditions 110|) ; let us define D(z\,z) by 

c f z dz 1 

D(z 1 ,z) = —(l + zi)D(z 1 )D(z) 



Hn 



D 2 (z')g(z') 



, (11) 
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subject to the conditions 

D(zi,z)\ z = Zl =0, 

d , . , c 1 + zi 

— L> Zi,z) = Slgn {z-zi)— — . 

dz H g(zi) 



(12) 



We see that D(z±, z) satisfies Eq. if the function g(z) 
is a solution of the equation 



lng(z) 



dz 
so that 



d-z 
dv 2 



dv J 



r d2z ', 

g(z) = g exp / -^^dz' 
J (it) 



(13) 



(14) 



where go is an arbitrary constant of integration. It is easy 
to relate the function g(z) to the Hubble function H(z). 
In fact, using 



dz 

~dt 

and 

dz 



= -(l + z)H(z), 



,H(z) 



d~ V =( 1 + Z Y H ' 
it follows that Eq. I|13[) can be rewritten as 

.H(zY 



^-\ng{z) = ^~ 

az dz 



ln(l + z) 



H 



so, 



g(z) = (l + z) 



a g(f) 



(15) 



(16) 



(17) 



(18) 



where we have imposed the initial condition g(0) = 1. Let 
us remind that one of the Friedman cosmological equations 
can be rewritten as 



h(z) = H y/n m (i + zf + n fc (i + z y + n A , 



(19) 



where Ho is the present value of the Hubble constant, 

8irGg m . , 

= is the matter density parameter, ilk — 

3Ho 

kc 2 , , n Ac 2 . , 

is the curvature parameter and s2a = -7-^- is the 



R0H0 3 Ho 

cosmological constant parameter. The omega parameters 
satisfy a simple algebraic relation 



fl m + flk + = 1- 



(20) 



The initial conditions l|13|) expressed in terms of H(z) be- 
come 



D{z\,z)\ i 



0, 



d c 

—D(z u z)\ z=Zl ^sign(z-zi)— . 

dz {l + Zi)H(zi) 



(21) 
(22) 



2.3. The angular diameter distance 

To apply the angular diameter distance to the re- 
alistic universe it is necessary to take into account 
local inhomogeneities in the distribution of matter. 
Unfortunately, so far an acceptable averaging procedure 
for smoothing out local inhomogeneities has not been 
developed (Krasinski 1997). Following Dyer & Roeder 
(Dy er fc Roeder 1972| ), we introduce a phenomenological 



parameter a 



1 



Bclv 



<e> 



called the dumpiness pa- 



rameter, which is related to the amount of matter in 
clumps relative to the amount of matter distributed uni- 
formly l|Dashveski fc Zeldovich 19651 ITomita et a!.1999|l . 
Therefore, Eq. 0, in the FRW case, can be rewritten in 
the form 

dz\ 2 d 2 D ( d 2 z\ dD 3 „ , . r 
-I — + f — I — + -att m {l + zfD = Q. (23) 



dv J dz 2 



dv 2 I dz 



Let us note however that this equation does not fully de- 
scribe the influence of non homogeneities in the distribu- 
tion of matter on the propagation of light. It takes into 
account only the fact that in locally non homogeneous uni- 
verse a light beam encounters less matter than in the FRW 
model with the same average matter density. Equation 
(|23|l does not include the effects of shear, which in gen- 
eral is non zero in a locally non homogeneous universe 
(see, for example SEF Chapter 4 and 11). The effects 
of shear have been estimated with the help of numeri- 
cal simulations (see, for example Schnei der fc Weiss~1 988 
ILinder 19981 |Jaroszynski 1991[TWatanabe fc Sasaki 199f)j> 
but satisfactory analytic treatment is still lacking. 

It is customary to measure cosmological distances in 
units of c/ Ho, therefore we introduce the dimensionless 
angular diameter distance r = DH /c. Using Eq. I|16fl 
and Eq. H19|) . we finally obtain 

(1 + z) [n m (i + zf + n fc (i + z) 2 + n A ] + (24) 



7 \ dr 

-n m (i + zf + 3Ji fc (i + zf + 2n A — 

2 ) dz 



+-aQ m (l + zfr = 0, 



with initial conditions 

K*)U=o = o, 



dr(z) 
dz 



z=0 — 1- 



(25) 



(26) 



Eq. H24f) can be cast into a different form by using u = 
1/(1 + z) as an independent variable instead of z; we get 



d 2 r 3 dr 
u 2 (tt m + Q k u + Hau 3 )—^ - u(-Q m + fl k u) — 



du 2 



dv 



--aCt m r = 0. 



(27) 



On the other hand, using the cosmic time t as an inde- 
pendent variable, Eq. (|23|l assumes the form 



dt z dt 



(28) 
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This equation was for the first time introduced by 
Dashevski & Zeldovich HDashveski fc Zeldovich 1965)) (see 
also |Dashveski fc Slysh 1966| ) , and more recently Kayser 
et al. ( |Kayser et a!.1997| ) and IjKantowski 1998|l have used 
it to derive an equation similar to Eq. In Eq. 

the dumpiness parameter a is usually considered as a 
constant. However, in IjDashveski fc Zeldovich 1 965 and 
|Kayser et al.l997| , a is allowed to vary with time, but 
only the case a = constant is really considered. For a dis- 
cussion of the general case when a depends on z see, for 
example, the paper by Linder IjLinder 1988|l . 
To give an example of our procedure of constructing the 
angular diameter distance between a fictitious observer at 
zi and a source at z 2 , let us consider the simple case of 
a = constant, Q\ = 0, and f)j. = 0. In this case it is easy 
to integrate the equation (22), we obtain 



r(z) 



{I + zf - {I + z)- 
2/3(1 + z) 5 / 4 



(29) 



where (3 = \\/2h — 24a, is the general solution 
of Eq. (|25|) with the imposed initial conditions. 
This solution can be found, for example, in SEF 
( |Schneider, Ehlers fc Falco 1992| ), (see Eq. (4.56) there). 
The function g(z) can be easily computed in the flat A = 
Friedman-Robertson- Walker cosmological model; in fact, 
in this case we have 



H(z) = H (l + z)Vz + l, 

so that, according to the Eq. H18|) . we have 

g{z) = (1 + zfs/z + 1. 

Thus using (|ll(l . we get the familiar solution 



D{z U z s ) = - 



found in SEF. 



(i + z 2 f-i (l + z^-i 

(l + z^+i (1 + 31)0+4 



(30) 
(31) 

(32) 



3. Exact solutions 

In the general form the Eq. I|24|l is very compli- 
cated. From the mathematical point of view this 
equation is of Fuchsian type with four regular sin- 
gular points and a regular singular point at infin- 
ity Qlnce 1964|l . General properties of this equa- 
tion have been extensively studied by Kantowski 
UKantowski 19981 |Kantowski, Kao fc Thomas 2000| 
IKantowski fc Thomas 2001|> . 

Following Kantowski let us cast the Eq. I|24fl into a 
different form replacing 1 + z — x and introducing 
r ( z ) — T+i'; we obtain 



(fir. 



Qi.x 2 



-(a - l)ft m x 



„ N d 2 h /3_ 



Q k )h = 0. 



lh I x— 

ax 



(33) 



When Qk ^ by rescaling x, x = jy^x this equation can 
be turned into the Heun equation 



, 2 ^A^m.^h .3, ,,~dh 

{x +x + ^f ) d^ + { r +1)x Tx 



— [a 

l 2 y 



l)x - 11 ft. = 0. 



(34) 



The angular diameter distance can be expressed in terms 
of basic solutions of the Heun equation as 



r(z) 



ft(^(l 



(35) 



provided that h(%£) = and h'(jQ-) = The generic 
solution for the angular diameter distance expressed in 
terms of the Heun functions is very complicated, in an 
explicit form it is given by Kantowski IjKantowski 1998)l . 
When ilfc = 0, dividing Eq. (|33|l by f2 m , we obtain 



, o £l^,d 2 h 3 9 dh 3. , 

( * + sn ) dx^ + rd-x + 2 ia - 1)xh = °- 



By changing variables to rj{x) 



1 



(36) 



and intro- 



ducing new independent variable P(x) = x 3 ^ 4 h(x) this 
equation can be transformed into the associated Legendre 
equation 

d 2 P 



(!-„')- 



2i) 



dP 



1 



dr/ 2 

In this case (ilfe 
given by 



(5- 



25 - 24a 



)P = 0. 



(37) 



dr\ 36 v (1 — 7] 2 ) 

0) the angular diameter distance is 



T^x (38) 




where S = \/25 — 24a, and P£(z) denotes the associated 
Legendre function of the first kind. In section 5 we will 
propose a simple approximate analytic solution of the gen- 
eral equation for the angular diameter distance which is 
reproducing reasonably well the exact numerical solution 
in the range of redshifts (0, 10) far exceeding the range of 
redshifts of observed supernovae of type la. 

4. Angular diameter distances in the gravitational 
lensing theory. 

The angular diameter distances, or their combinations, 
appear in the main equations and in the most important 
observational quantities of the gravitational lensing (GL) 
theory. This fact, together with the relation between the 
angular diameter distance (Da) and the luminosity dis- 
tance (Dl) 



D L = {l + z) 2 D A , 



(39) 
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makes the study of the equation for the angular diameter 
distance still more important, also from the point of view 
of better interpretation of observational data. In this sec- 
tion we describe very briefly the role of angular diameter 
distance in the gravitational lensing theory, mentioning 
only some effects in which it plays an important role. 
Let us begin with the general expression for the time delay 
between different light rays reaching the observer 



cAt = (1 



Zd) 



DdDs 
D, 



ds 



(8 - py - } + constant, (40) 



where Dd, D s and Dd s are correspondingly the angular 
diameter distances to the deflector, to the source and be- 
tween the deflector and the source. The term (0 — (3) in 
Eq. H40f) represents the geometrical time delay and the 
other term is connected with the non homogeneous dis- 
tribution of matter ( |Schneider, Ehlers fc Falco 1992| ). The 
first term in Eq. I|4(J|) is proportional to an important com- 
bination of angular diameter distances namely, to — - — - . 

Dds 

This combination of angular diameter distances was in- 
troduced for the first time by Refsdal (|Refsdal 1966J) . and 
was used to obtain an estimate of Hq from the time delay 
measurements of multiply imaged quasars. It is possible 
to rewrite it in the following way 



(1 + z d ) 



Dd, 



-no 



(41) 



where the function \ is given by 
f°° dz 

x(z,n m ,tt A ,Q k ,a) = / 2 -, , 
J z r 9(z) 



Ho 



r 2 (z)(l + z) 2 H(z)dz' 



(42) 



So, x is connected with the general solution of Eq. 
and it directly appears in the expression for time delay. 
Let us now consider the cosmological lens equation. The 
metric describing FRW cosmological models is confor- 
mally flat. In this case the simplest way to derive the lens 
equation is to use the Fermat principle, so we have 



dAt 
or 



0. 



2R 



dip 



P = - + Zd){x(zd) - 



cH 



dd 



(43) 



(44) 



where R„ is the Schwarzschild radius of the deflector. By 



denoting^ = DdO, rj = D s /3 7 and a(£) 
transform Ea . (|44J) into 

L>d 



2Rg 

'Dd'oe' 



(45) 



which is formally identical with the lens equation for 
z << 1. As is apparent, in the lens equation (|45|l another 
dimensionless combination of angular diameter distances, 



D s /Dd appears, besides the angular diameter distance it- 
self, Dds- In other words, in the equations describing grav- 
itational lensing we find quantities, which can be written 
in terms of the angular diameter distance D(z) and the 
X{z) function. In the general case D{z) and x( z ) can b e 
evaluated only numerically. In the next section we propose 
a simple analytical approximation for both functions. We 
hope that these approximate forms can be useful, for ex- 
ample, in big numerical codes used to derive basic cosmo- 
logical parameters from observational data and to study 
for instance, statistical lensing, weak lensing, microlensing 
of QSO, etc. 

5. The approximate expression for the angular 
diameter distance 

In the generic case the equation l|2^|) does 
not have analytical solutions (Kantowski 1998 
|Kantowski, Kao fc Thomas 2000| ). From the mathe- 
matical point of view this equation is of the Fuchsian 
type fince 1964fl with four regular singular points and a 
regular singular point at infinity. The solutions near each 
of the singular points, including the point at infinity, 
are given by the Riemann "P-symbol and in general the 
solutions can be expanded in a series of hypergeometric 
functions IjTricomi 1961J1 or expressed in terms of the 
Heun functions IjKantowski 1998)l . 

In practice, in the general case, the equation (|24|l with 
appropriate initial conditions is solved numerically. In 
Fig.l we show a numerical solution of the D-R equation 
for a = 0.8, n A = 0.7, fl k = and Cl m = 0.3. Using 




Fig. 1. The dimensionless angular diameter distance as a 
function of the redshift in a flat cosmological model with 
n m = 0.3, n A = 0.7 and a = 0.8. 



Mathematica we discovered that there is a simple func- 
tion r(z) which quite accurately reproduces the exact 
numerical solutions of the equation 124(1 for z up to 10, it 
has the form 



r(z) 



\Jd lZ 2 + (l + d 2 z + d 3 z 2 Y 



(46) 
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where d\ , d 2 and are constants which depend on the pa- 
rameters that specify the considered cosmological model. 
First of all please note that the function i|46|) automati- 
cally satisfies the imposed initial conditions, so r(0) = 
dv 

and -^-(0) = 1. To express the constants d\, d 2 and c?3 

dz 

in terms of f2A> ^m, ^fc and a we inserted the proposed 
form of r(z) into the equation (|24[) and required that it 
be satisfied to the highest order. In this way we obtained 
that: 



d 3 



(2 fe + 3fi ro ) 2 7Q k + (13 + a) n r , 



14 n k + (28 + 5 a) fi„ 
6 Q k + 7 Sl m + 4 Q A 

2 ^^/j; | 3 ^if-yi 



1 



l -((2n k + m m f - A (50 fe + 9^ 
.6 V 



16 

8 (14^ + (28 + 5a) Q m ) 

6fi fe + 7n m + 4n A 



(47) 

(48) 
(49) 



In Fig. 2 we show the approximate solution of the equation 
(EU) for a = 0.8, n A = 0.7, = 0, and Q m = 0.3 and for 
comparison we also plot the exact solution. 



the exact numerical solution for xi z ) and for comparison 
we plot the approximate solution with the best fit value 
of 7 = 6. 



10 



2 

z 



Fig. 3. The exact dimensionless function \ (solid line) for 
a model with a = 0.8, fl\ = 0.7, fl k = 0, and fi m = 0.3. 
With dashed line we plot the analytic approximation with 
the best fit value of the parameter 7 = 6. 



. 4 



0.2 



. 1 



Fig. 2. The exact dimensionless angular diameter distance 
r is plotted (solid line) and compared with the analytic 
approximation (dashed line) in a model with a = 0.8, 
£1a = 0.7, Vt k = and Q m = 0.3. In the considered range 
of redshifts the relative error is less than 10%. 



Using Mathcmatica we have also found a useful ap- 
proximate analytical form for the x function 



xO) 



1 



1 



(50) 



where 7 is a constant. Unfortunately we have not been 
able to analytically relate 7 to a and other cosmological 
parameters and the appropriate value of 7 should be ob- 
tained by the standard fitting procedure. In Fig. 3 we show 



6. Conclusions 

In this paper we discuss the angular diameter distance 
in the Friedman-Robertson- Walker cosmological models 
and consider the case when the cosmological constant 
and the curvature of space could be different from zero. 
The effects of local non homogeneous distribution of mat- 
ter are described by a phenomenological parameter a, 
consistently with the so called empty-beam approxima- 
tion. Unfortunately, at the moment there are no gen- 
erally accepted models that describe the distribution of 
baryonic and dark matter and therefore the influence of 
inhomogeneities of matter distribution can be included 
only at this approximate level. In the generic case the 
equation l|24J) is of a Fuchsian type, with four regular 
singular points and one regular singular point at infin- 
ity. The general solution of this type of ordinary dif- 
ferential equation is given in terms of the Heun func- 
tions l|Kantowski 19981 IDemianski et al. 2000| l. 

However the exact solution is so complicated that 
it is useless in practical applications IjKantowski 1998 
|Kantowski, Kao fc Thomas 2000| ) . Therefore we have pro- 
posed an approximate analytic solution, simple enough to 
be used in many applications and at the same time it is 
sufficiently accurate, at least in the interesting range of 
redshifts (0 < z < 10). In Fig. 2 we compare the exact 
numerical solution of the equation for the angular diame- 
ter distance with the approximate one. The approximate 
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solution reproduces the exact curve quite well and the rel- 
ative error does not exceed 10%. 

Following SEF, we have found the function x which ap- 
pears in the expression for time delay as well as in the lens 
equation, and which naturally appears in the expression 
for angular diameter distance between two arbitrary ob- 
jects at redshifts z\ and z 2 ( see Eq. dJ). 
We have also proposed an approximate analytical form of 
the function \ which depends only on one parameter 7 but 
unfortunately 7 has to be fixed by a standard fitting pro- 
cedure (see Fig. 3). Our approximations have been already 
applied in complicated codes used to study the statistical 
lensing UPerrotta et al.2001|) . 

Finally we would like to stress that from our analysis it 
follows that variations in the angular diameter distance 
caused by the presence of cosmological constant are quite 
similar to the effects of a non homogeneous distribution 
of matter described here by the dumpiness parameter a. 
In Fig. 4 we plot the angular diameter distance for two 
models, one with a homogeneous distribution of matter 
(a = 1) and A / and another with a non homoge- 
neous distribution of matter (a = 0.7) and A = 0. We see 
that an inhomogeneous distribution of matter can mimic 
the effect of a non zero cosmological constant. This is 
an important observation in view of the recent conclu- 
sions based on observations of high redshift type la su- 
pernovae that the cosmological constant is different form 
zero (|Perlmutter 1997MRiess 2000)l . 
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Fig. 4. The angular diameter distance for two models, one 
with a homogeneous distribution of matter [a = 1) and 
A /: (solid line) and another with a non homogeneous 
distribution of matter (a = 0.7) and A = (dashed line). 
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8. Appendix 

As was already noted the equation for the angular di- 
ameter distance in a locally non homogeneous Friedman- 
Robertson- Walker universe with non zero cosmological 
constant can be reduced to the Heun equation (see Eq. 
28) (Heun, 1889). The Heun differential equation general- 
izes the Gauss hypergeometric equation, it has one more 
finite regular singular point (see Ince, 1964). This equation 
often appears in physical problems in particular in study- 
ing of diffusion, wave propagation, heat and mass transfer, 
and magnetohydrodynamics (see Ronveaux, 1995). In the 
general form the Heun equation can be written as 

d 2 H ( a,\ (X2 83 \ dH 

dy 



d y \y-yi y-ni 

(a-4y - q) 



y - V3 

H = 0, 



(51) 



{y - yi){y - yi){y - yz)' 

where a\, a-i, 03, 04, q, yi, 1/2 and 2/3 are constants, or in 
a canonical form as ( |Bateman <k Erdelyi 1955| ) 

/ \ / \ d 2 H 
x{x — l)(x — a)- 



dx 2 



{(a + /3+ l)x 2 - 
{a/3x - q)H = 0, 



dH 



+ 13 + 1 + 0(7 + S) - 6] x + 07} — + 

X (52) 



where a, (3, 7, <5, a, and q are constants. 

When 7 is not an integer, the general solution of the 
Heun equation can be written as: 

H(x) = C 1 F(a,q,a,0,'Y,5,x)+ (53) 
CM^Fia, q u a - 7 + l,/3 - 7 + 1, 2 - 7, 5,x), 

where C\ and C2 are constants and qi = q + (a — 7 + 
l)(/3-7+ 1) -a/3 + 6(1- 1). 

Following our success in finding an approximate so- 
lution of the equation for angular diameter distance we 
would like to present an approximate solution of the Heun 
equation. It has the form: 



H(x) = B- 



xo \xq 



(54) 




1 



1 



where the parameter xq represents the point where the 
initial conditions are specified and di, d%, d%, B are con- 
stants. For the Heun equation we adopt the following ini- 
tial conditions 1 : 

H[x ] = 0, (55) 

1 To get the general initial conditions i.e. H(xo) = A it is 
enough to add to the 155H a constant H + A. 
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dH 



B 



(56) 



Please note that the form 1|55|) automatically satisfies the 
imposed initial conditions. 

To express the constants d\, d,2 and c?3 in terms of a, 0, 
7, S, a, and q, it is necessary to insert the proposed form 
of H (x) into the canonical Heun equation 1|52|) and require 
that it be satisfied to the highest order. The constants d\, 
c^2 and g?3 obtained in this way are quite complicated. It 
will take about 3 pages to give them here. They are shown 
in an explicit form at |http: / / people. n a.infh.it~ester| 




Fig. 5. Comparison of the exact (solid line) with the ap- 
proximate solution (dashed line) of the Heun equation. 
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Fig. 6. Relative error of the approximate solution of the 
Heun equation as a function of x. 



In Fig. 5 we compare an exact numerical solution of the 
Heun equation with the approximate solution for the same 
values of a, 0, 7, 6, a and q. 
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